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1). 

Abstract: We study holographic models describing an RG flow between two fixed points 
driven by a relevant scalar operator. We show how to introduce a spurion field to restore 
Weyl invariance and compute the anomalous contribution to the generating functional in 
VjD '■ even dimensional theories. We find that the coefficient of the anomalous term is pro- 

portional to the difference of the conformal anomalies of the UV and IR fixed points, as 
expected from anomaly matching arguments in field theory. For any even dimensions the 
J"""- 1 coefficient is positive as implied by the holographic a-theorem. For flows corresponding to 



spontaneous breaking of conformal invariance, we also compute the two-point functions of 
the energy-momentum tensor and the scalar operator and identify the dilaton mode. We 
find that the simplest models with one scalar field have IR divergence which is at odds 
with unitarity in a local field theory. 
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1. Introduction and summary 

Recently Komargodski and Schwimmer have proved the four dimensional "a-theorem" [1,2]. 
This theorem, conjectured originally by Cardy [3], states that for any Renormalization 
Group (RG) flow between an UV and an IR fixed points the following inequality holds 

auv > am (1-1) 
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where auv an d am are the UV and IR conformal anomalies related to the expectation 
value of the energy-momentum tensor as follows 

<r;> = cW 2 - aE A , (1.2) 

where W 2 is the Weyl tensor squared and £4 is the Euler density. The a-theorem is the 
four dimensional relative of the well known two dimensional Zamolodchikov's c-theorem [4]. 
The proof of the a-theorem of [1] follows from the analyticity and unitarity of the forward 
4-point scattering amplitude of a scalar field t(x). The latter corresponds to the dilaton, 
the Goldstone boson associated with the spontaneously broken conformal invariance or to 
the conformal compensator (spurion) for RG flows driven by relevant operators. Another 
interesting result based on the ideas behind the a-theorem is a proof that perturbative 
theories always flow to CFTs [5] , although it seems that some aspects remain to be under- 
stood [6]. 

Since there is a monotonically decreasing function associated with RG flows in two and 
four dimensions, a natural question to raise is whether this theorem can be generalized to 
field theories at any even dimensions. It turns out that the generalization of the arguments 
of [1] to six dimensional field theories is not straightforward as was revealed in [7]. An 
alternative to proving a generalized theorem directly for d dimensional field theory is to 
do it using a dual d + 1 dimensional holographic gravitational bulk system. This is one of 
the goals of the present paper. 

RG flows have been intensively studied using holography. In most of the cases (e.g. 
[8-15]), the RG flows are driven by a relevant deformation of the Lagrangian. RG flows 
associated with spontaneous breaking of conformal symmetry were studied in, for in- 
stance, [16-18]. Based on holography, a-theorems in arbitrary number of dimensions where 
proven using various different approaches. By studying the renormalization group flow 
along null geodesic congruences the authors of [19] have proven a holographic version of an 
"a-theorem" . An "a-function" that is positive and monotonic if a weak energy condition 
holds in the bulk gravity theory was constructed in [10]. For even-dimensional boundaries, 
the a-function was shown to coincide with the trace anomaly coefficients of the correspond- 
ing boundary field theory in limits where conformal invariance is recovered. In [20] it was 
demonstrated that unitary higher curvature gravity theories in both odd and even dimen- 
sions admit a-theorems in terms of their corresponding entanglement entropy. In cases 
where the dual CFT is even-dimensional it was shown that it is the coefficient of the "a 
anomaly" that flows. 

The RG flows we study have a holographic dual description as a solitonic solution 
interpolating between two AdS spaces of dynamical Einstein gravity coupled to a single 
scalar field in (d + 1) dimensions, which is generated by an appropriate scalar potential. 
We determine the solution of the scalar field and of the metric that admits a d-dimensional 
Poincare isometry. The scalar and the metric both depend only on the (d+l) th dimensional 
(radial) coordinate. We then introduce the Fefferman-Graham coordinates to express the 
general fluctuations of the metric components and of the scalar field in the AdS regions in 
terms of a power expansion in the radial coordinate which is also an expansion in powers 
of the d dimensional space-time derivatives. This expansion is valid up to some value of 
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the radial coordinate, so we need to introduce a cutoff in the radial direction that acts as 
an IR regulator. We then introduce the spurion via a generalization of the Penrose-Brown- 
Henneaux (PBH) diffeomorphism [21,22], whose effect is to perform a Weyl transformation 
on the boundary values of the metric and the scalar. The same kind of transformation (in 
pure AdS) has been used recently to compute the effective action of the dilaton from the 
action of probe branes [7] . A generalization of the transformation to theories with scalars 
was also used in [23] to compute contributions of relevant operators to the entanglement 
entropy. Let us stress here that in our case there is no analog to the scalar fields on the 
probe branes and the parameter of the PBH transformation is not a dynamical mode. 1 

The PBH transformation admits the same kind of expansion as the metric and scalar 
fluctuations, in such a way that terms in the expansion transform covariantly. There 
is an exception that appears at dth order in derivatives, whose effect is to introduce a 
term related to the Weyl anomaly in the generating functional computed using the usual 
holographic prescription. The term takes the following form 

•d„ I ~t,^(~\ ( aUV ur 



d d x^-gr(x) (A u d v ~A d H , (1.3) 



where t(x) is the spurion field and g is the determinant of a background metric, dressed 
with the spurion to make a Weyl-invariant combination. Ad is the conformal anomaly in 
even dimensions, at the UV or IR fixed points. In general, the anomaly term takes the 
form [24] 

A d = -aE d + J2 c iIi, (1.4) 

where E d is the Euler density and Ij are Weyl-invariant polynomials of the curvature and 
its derivatives with d derivatives of the metric, examples in d = 2 and d = 4 are given 
by (|3.5| ),( ^6|) . When the action of the gravitational dual is Einstein gravity plus matter 
fields, as in our case, the value of the coefficients q is proportional to a with fixed factors. 
In more general cases than the ones we consider the coefficients q are independent. Our 
computation confirms that there is a term in the generating functional of the form 

-(a uv - a IR ) J d d x a/^t(x) E d . (1.5) 

In deriving this result we have used the fact that 

a uv -a IR ^\ Uuv 11 - 1 - Lm*- 1 ) , (1.6) 

with Ljjv an d Lir denoting the radii of the two asymptotically AdS spaces and k is the 
d dimensional Newton constant. Since according to [10] Ljjv > Lm, 2 it indeed follows 
that the coefficient of the anomalous term auv ~ a iR > is positive. This result is in full 
accordance with the "a-theorem" of [1]. We would like to emphasize that whereas in field 



1 However, it could be promoted to a dynamical mode by changing the kind of boundary conditions of 
the metric, we discuss this possibility in Sec. 



4.e 



For this to be true, the enery-momentum tensor of the bulk matter fields should satisfy the null energy 
condition Tmnu m u N > 0, with qmnu u u n = 0. 
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theory such a statement has been proven in two and four dimensions, our result holds at 
any even dimension. 

In addition to the study of the holographic "a-theorem", we analyze in detail the 
manifestations of the Goldstone theorem related to the spontaneous breaking of conformal 
symmetry by computing the one and two point functions of the operator dual of the scalar 
and the energy momentum tensor. To implement this, we consider a suitable potential for 
the scalar field which admits the appropriate asymptotic behavior in the UV correspond- 
ing to a vacuum expectation value (vev) of a dual operator O of the CFT\jy. We first 
expand Einstein's equations to linear order in the fluctuations, fixing a particular (radial) 
gauge. We then analyze the tensor and scalar fluctuations at low momenta away from the 
horizon and compare with the finite momentum solution close to the horizon. By carefully 
matching the two asymptotic behaviors we can determine the full solution to leading order 
in momentum. Evaluating the on-shell action on these solutions we can compute corre- 
lation functions using the standard holographic prescription. First, we confirm that the 
operator dual to the scalar field has a non-zero expectation value (0), while the trace of 
the energy-momentum tensor vanishes. The two-point function of the energy- momentum 
tensor reads 

L IR d 1 



T>»{-q)T#{q)) oc > H U^(^) d log^y^) 2 , (1.7) 

where is the four momentum and n^'"' 3 is a kinematic factor depending on the mo- 
mentum and the metric which is transverse and traceless. The explicit expression is given 



in (4.g). The two-point function involving the energy- momentum tensor and the scalar 



(T^(-q)0(q)) = |^M-Lp^ - (O) ^ . (1.8) 



operator is given by 



where PJ^ U is a kinematic factor transverse to the momentum (see ( [4.8D ). This result is 
in agreement with the field theory expectation based on Goldstone's theorem. For the 
two-point function of the scalar operator we find 



K 2 LlR d+2-2A IR 



(0(-q)0(q)) OC ^ (O) 2 (v^) d -^n. (1.9) 

This result exhibits a long range interaction, which is consistent with Goldstone's theorem 
(a more detailed explanation is given in Sec. |4.5| ), but the zero-momentum singularity is not 
a pole, as one would expect to find for a free dilaton mode, and in fact the singularity leads 
to an IR divergence that is at odds with unitarity. Recall that naively one would expect to 
find in the IR a CFTm x free dilaton. We comment more on this in the Discussion section. 
The absence of a pole is in contrast with the probe brane analysis [7], where the dilaton is 
a mode which essentially decoupled from the remaining massless degrees of freedom. 

The paper is organized as follows. In Section ^ the basic setup of the holographic RG 
flows which is analyzed in this paper is established. Section [3| is devoted to the manifestation 
of the conformal anomaly matching in holography. We first briefly review the picture of the 
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boundary field theory. A compensator field (spurion) is introduced to maintain conformal 
invariance of the field theory coupled to an external curved space-time. The anomalous 
terms of the generating functional are written down and the corresponding "a-theorem" [1] 
is stated. Back to the holographic dual, we show that the generating functional of the field 
theory includes an anomaly term whose coefficient is given by the difference between the 
a-anomaly of the UV and of the IR, in full accordance with the "a-theorem of [1]. The 
holographic spontaneous breaking of conformal invariance is discussed in Section ||. We 
determine the tensor and scalar fluctuations at low momenta and extract the holographic 
correlators (O) , (OO) , (T^O) and (T^T a/3 ). In Section [D| we promote the spurion field 
to a dynamical mode by introducing a new set of boundary conditions for the metric. We 
conclude and discuss several open questions in section || We end this paper with four 
appendices. In the first we elaborate on the equations of motions of the fluctuations. In 
appendix [B] we elaborate on the matching procedure and explicitly show the existence 
of an overlapping region between the boundary and the near-horizon region. Appendix 
[C] is devoted to an application of our matching procedure to the case of the Coulomb 
branch of the N = 4 SYM theory that was previously analyzed using a different procedure 
in [15,25-27]. In appendix [i^ we present a family of toy models that describe an RG flow 
between two fixed points. 



2. Holographic RG flows 

A holographic dual of a theory with an RG flow between two fixed points will be a geom- 
etry that interpolates between two anti-de Sitter spaces with different radii. In Gaussian 
coordinates, the metric will be 

dsl v = dr 2 + e 2r/Luv i] llu dx fl dx' / , r -)• +oo, 

dsj R = dr 2 + e^^'n^^dx^dx", r -> -co, (2.1) 

where, imposing some physical conditions, the radius close to the boundary r — > oo is 
larger than the radius close to the horizon r — > — oo: Ljjy > L/_r. The radial coordinate r 
maps to the RG scale of the dual theory, while the coordinates x^ map to the space-time 
where the dual theory lives. 

Asymptotically the geometry has the following isometry 

r^r + A, x 1 * -»• e~ A /V*, (2.2) 

where A is a constant and L = Ljjv oi L = Ljr depending on which limit we take. A 
shift in the r coordinate is then equivalent to a dilatation transformation of the space-time 
coordinates. 

There are examples in supergravity of geometries that interpolate between two different 
AdS spaces (e.g. [8-14]), but they describe RG flows driven by a relevant deformation of the 
Lagrangian and not spontaneous breaking of conformal symmetry. Since we are interested 
mainly in the latter, we will construct our own holographic toy models with the desired 
properties, keeping the discussion as general as possible, so the results could easily be 



- 5 - 



extended to particular examples in supergravity in case some are eventually found. Our 
construction is inspired by similar models that have been studied in the past [16-18]. 

It is relatively easy to construct holographic models dual of a flow between fixed points, 
taking as inspiration some of the known supergravity flows. The simplest case is a scalar 
field with some potential V coupled to Einstein gravity in d + 1 dimensions. 



Seh = ^ f d d+l x^ (~R - X -d M <\>d M 4> ~ V{<f>) ) . (2.3) 



The requirements on the potential are 

a) There are at least two critical points where the potential is negative, one is a maximum 
and the other is a minimum. 

b) There is a classical solution that interpolates between the two critical points. 

The general case is still too complicated, but in certain cases it is possible to find a system 
of first order equations whose solutions are also solutions of the full second order system. 
That is the case when the potential is written in terms of a superpotential W s 



2 



(2.4) 



One introduces an ansatz for the metric with Poincare invariance along the space-time 
directions 

ds 2 = dr 2 + e 2A ^r,^dx»dx v , (2.5) 

and the first order equations that give solutions for the Einstein plus scalar equations of 
motion are 

<t>' = -dW, A' = -^, (2.6) 

where primes denote radial derivatives. Critical points of the superpotential are also critical 
points of the potential, so in order to have a flow we should demand that the superpotential 
has two critical points. Without loss of generality we can assume that one of the critical 
points is at 4> = 0, 



W^^l + gTLf, v ^_ d(d-l) + A uv (A uv -d) (j)2 
Ljjv 2Luv 2Ljjv 2Luv 



Here Ajjv is the dimension of the operator dual to the field (p. We will assume in the 
following that the dual operator is relevant, so Ayy < d. In this case the potential has a 
maximum. 



3 Locally in field space the potential can always be written in terms of a superpotential, by solving 



the equation (2.4) understood as a differential equation. For each solution one has a different RG flow 
geometry, that could be dual to a theory with explicit or spontaneous breaking of conformal invariance. 
Here we assume a prior knowledge of the superpotential, which is assumed to be valid in all configuration 
space, but the analysis could be extended to situations where this is not the case. We thank Ioannis 
Papadimitriou for pointing this out. 
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4> <P 

Figure 1: An holographic RG flow interpolates between the minimum and the maximum of the 
superpotential on the left side, or equivalent ly between the maximum and the minimum of the 
potential on the right side. 



The second critical point will be at some value <f) m , such that 

LlR ^IR ILlR I 



We will demand that Ajr > d, so the dual operator becomes irrelevant in the IR fixed 
point and the potential has a minimum. Solutions can be simply obtained by using the 
field <j) as radial coordinate, 

A family of simple superpotentials satisfying the requirements is given in Appendix [D|. In 
Figure [| we have plotted an example of a superpotential with the required properties and 
its corresponding potential. 

The construction can be generalized to several scalars, that may have a non-trivial 
metric in configuration space. The basics will be the same, the flow between two fixed 
points will be dual to a solution to the equations interpolating between two critical points 
of the superpotential with the right properties. In this case one has to worry about the 
directions transverse to the flow in configuration space, and make sure that there are no 
unstable modes. 



3. Conformal anomaly matching in holography 

Our purpose in this section is to show explicitly how the anomaly matching arguments 
presented in [1,2] for the generating functional are realized in gauge/gravity models. 

3.1 Field theory preliminaries 

Consider first a conformally invariant theory. In the presence of a background metric 
conformal invariance is explicitly broken. However, one can introduce a compensator field 

4 See also [28] for the original arguments regarding anomaly matching in theories with spontaneously 
broken conformal invariance. 
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such that the theory is invariant under an extended symmetry. The compensator field, or 
spurion, will enter as a Weyl factor of the background metric, as was discussed in [1,2], 



9nv -> 9fiu = <? T 9iLv (3.1) 

The extended symmetry involves both a Weyl rescaling of the metric and a shift of the 
spurion: 

9iiv ->• e 2a 9ij,v, r^T-a. (3.2) 

Correlation functions of the energy-momentum tensor in the CFT can be obtained by 
deriving a generating functional (generally non-local), T, with respect to the background 
metric. The generating functional should have the same symmetries as the original theory, 
so it should be a diffeomorphism-invariant functional of the dressed metric g^ u . In the 
absence of anomalies this will be the end of the story. However, when the number of 
dimensions d is even, a CFT in curved space has in general a conformal anomaly. The 
generating functional should reproduce the anomaly, which implies that it is not completely 
invariant under the transformation (13.21) but it transforms as 



<5 CT r| T=0 = / d d x^aA d , (3.3) 



where Ad is a polynomial of the curvature and its derivatives with d derivatives of the 
metric. The anomaly in the generating functional coincides with the expectation value of 
the trace of the energy- momentum tensor. 



x-^W — 9iiv i c 



sr 



= Ad. (3.4) 

r=0 



For instance, in two-dimensions the anomaly polynomial is 

A 2 = -^R, (3.5) 

where c is the central charge of the CFT. In four and larger dimensions the number of 
independent coefficients is larger. Concretely, in four dimensions there are three possible 
independent contributions to the anomaly 

Ai = cW 2 - aE A + bOR, (3.6) 

where W 2 is the Weyl tensor squared and E± is the Euler density. The last term is not 
directly related to an anomaly, but it is a contact term that may appear in the energy- 
momentum tensor correlators and that can be eliminated through the addition of a local 
counterterm. The first and second terms are labelled by the coefficients c and a. The a- 
anomaly (known as type A) is the equivalent to the conformal anomaly in two dimensions, 
and it shares some properties with chiral anomalies: it is determined by a one-loop diagram 
with d/2 + 1 legs, there is no ultraviolet divergence involved and there is a set of descent 
equations related to the anomaly [24]. The c-anomaly (known as type B) on the other 
hand is related to UV divergences of correlation functions (see e.g. [28]). In dimensions 
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larger than four there is also a type A anomaly proportional to the Euler density. Type 
B anomalies also appear but their number increases with the number of dimensions, their 
classification is based on the construction of Weyl-invariant polynomials of the curvature, 
although in the presence of other sources there can be additional B anomalies. 

Now consider a situation where conformal symmetry has been broken, either sponta- 
neously or explicitly, and the theory flows to an IR fixed point. In the latter case it is 
still possible to maintain a symmetry of the form (3.2), by dressing dimensionful quantities 
with the spurion field. In other words, any mass scale m should be multiplied by a factor 
e~ T . In this way the generating functional in both cases takes essentially the same form. 
The generating functional will suffer in general from IR divergences due to the massless 
degrees of freedom of the IR CFT. One can avoid this by introducing an IR regulator /i 
and integrating out all the degrees of freedom above this scale, but keeping the degrees of 
freedom of the IR CFT. The resulting partition function will take then the form 

Z = e r ^^Z IRCFT [g^]. (3.7) 

With the IR regulator, F[fj,; g^] is a local functional of the dressed metric and couplings. 

From the anomaly matching arguments of [28] and [1,2], the variation of the total 
partition function is determined by the coefficients of the conformal anomaly in the UV 
theory, auv anci c uv- However, the IR CFT in general contributes to the anomaly with 
different coefficients, ajR and cir. Therefore, the generating functional r[/x;g MV ] must 
account for the difference. This implies that there is a term in the generating functional 
whose variation gives 

8 a T\ns g^} = J d d x^~ga {A u d v - A 1 /) . (3.8) 



For instance, in two dimensions 



T ^9,u] D - CUV 247T CIR J d 2 x^r (R+(dr) 2 ) . (3.9) 

The terms in four and six dimensions can be found in [1] and [7] respectively. In general, 
the coefficients of the anomalous terms coming from the variation of the functional should 
be proportional to the difference between the coefficients of the anomalies in the UV and 
IR CFTs. In two and four dimensions this allows to prove a c- or a- theorem [1,2], 

cuv > cir, auv > am. (3.10) 

Of course in two dimensions the c-theorem was proven long ago by Zamolodchikov using 
a different approach [4]. The new proof is based on the analytic properties of scattering 
amplitudes, unfortunately it seems that the extension of these arguments to six dimensions 
are not straightforward and a similar statement does not exist there [7]. 

This is in contrast with holographic models, where there are various existing proofs of 
a-theorems in arbitrary number dimensions [10,19,20]. We will show that the proof given 
in [10] for RG flow geometries implies that the coefficient of the anomalous term in the 
generating functional is positive. Therefore, there may still be an argument purely within 
field theory that fixes the sign of the coefficient in six dimensions, even though it would be 
more subtle than in two and four dimensions. 
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3.2 Spurion field for Weyl transformations 

We will now explain how to introduce a spurion field in a holographic dual. We start with 
the RG flow metric (2.5) and do the change of coordinates p = e~ 2A ^ 



do 2 1 
Ap 2 p 

where 



ds 2 = i 2 (p)^2 + zV^dx»dx v , (3.11) 



*(P) = (3-12) 

asymptotes Ljjy in the p — > limit and in the p — > oo limit. In these two limits the 
metric just becomes AdS in Fefferman-Graham coordinates. We can generalize this form 
of the metric to arbitrary solutions of the equations of motion 

ds 2 = £ 2 (p, x)^r + -g av{x, p)dx^dx u , (f> = <j>(x, p), (3.13) 
4p z p 

where, close to the boundary p — > the lapse factor becomes a constant £ — > Lyy. In 
this limit the metric and the scalar field will have an expansion in powers of the radial 
coordinate (for d even dimensions) [29] 

(x,p) = g$ (x) + g$ {x)p + ■ ■ ■ + g$ (x)p d ' 2 + h$ {x)p d / 2 log p 
+ ---+9 { l t Auv) (x)p d - Auv +■■■ 
<l>(x,p) = ( p(d-^uv)^ p (d-A uv )/2 + ... + ^Auv)^^^ + . . . ^_ 14 ^ 

For odd d dimensions, the logarithmic term proportional to hffi(x) is not present. As we 
will see this term is crucially related to the conformal anomaly. We are assuming that Aj/y 
is not an integer, otherwise there can be additional logarithmic terms in the expansion of 
the scalar field that will make the analysis more complicated. 

In this expansion all the coefficients are completely determined by the leading terms 
gffl and cf)( d ~ Auv \x), up to the terms gjfj(x) and (f)^ uv \x) for which the boundary 
conditions when p — > oo are also needed. The leading term in the expansion of the scalar 
field ^ d ~ A uv)( x ) jL o maps to a source for the dual operator in the field theory and breaks 
explicitly conformal invariance. The scalar solution still vanishes at the boundary, so the 
asymptotic AdS form of the metric is not affected. From the perspective of the holographic 
dual we are introducing a relevant deformation, so the UV physics will not be affected by 
it. 

(n) 

The coefficients g^J have been computed in pure gravity and in theories with scalar 
fields as well [26,29-31]. Each term has n derivatives, so the expansion above can also be 
seen as a derivative expansion. In particular, if the derivatives are small enough, it can be 
extended to values of p that are not asymptotically close to p = 0. This implies that for 
small enough derivatives 

\d a g$\, \d a ^\ «^ (3-15) 
Pir 

the expansion above will also be valid in a region of the geometry where the space asymp- 
totes AdS with radius Lm, that we will call the near-horizon region, even though a solution 
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with explicit dependence on the space-time coordinate will most likely continue to a dif- 
ferent geometry when p — > oo. 

Note that the solution in the near-horizon region is determined by the same leading 
term of the metric gffi , but that the scalar field will have a different expansion 

<f>{x,p) = (j) {d - AlR \x)p (d - AlR)/2 + ■■■ + ^ Air \x) P Air/2 + ■■■ (3.16) 

From the perspective of the near- horizon AdS, acts as a "source" although of 

course it is not independent of the boundary value 4>( d ~ Auv \ and the same is true for the 
subleading terms. For the explicit solution of the holographic RG flow, one sees that the 
scalar solution interpolates between the two scalings in a way determined by the superpo- 
tential. 

In the case where the geometry is the dual of a theory with spontaneous breaking of 
symmetry, the source term will vanish ^ d ~^uv) = q_ This is the case we have considered 
in the previous section, here we will keep the source term in order to make the discussion 
more general. 

We now proceed to introduce the spurion field. Note that in field theory the spurion 
was introduced by making a Weyl transformation of the metric and of dimensionful quan- 
tities. The question is what is the realization of this in the gravity dual. The answer for 
infinitesimal Weyl transformations in pure AdS is the so-called Penrose-Brown-Henneaux 
diffeomorphism [21,22] introduced a holographic context in [32]. We will generalize it to 
the RG flow geometries we introduced in the previous section and to non-infinitesimal 
transformations. The PBH transformation has also been used recently to construct the 
effective action of the dilaton from the action of probe branes embedded in AdS [7]. The 
difference with this approach is that the probe brane represents a breaking of conformal 
invariance that is suppressed in the large- N limit, so that to leading order the geometry 
that describes the UV fixed point is not affected and a decoupled dilaton shows up in the 
spectrum as a fluctuation of the brane. In the cases we study the breaking is not suppressed 
and the effect of the dilaton is seen indirectly in correlation functions. 

In pure AdS space the metric is (3.11) with a constant £ 2 (p) = L 2 . The infinitesimal 
PBH transformation takes the form 

p -> e~ 2r p ~ (1 - 2r)p, x^^-x^ + a* 1 , (3.17) 
a>1 = Y I" dp'g^{x,p')d v T(x), (3.18) 



o 



and it maps a solution with boundary metric gffl (x) to a new solution where the original 
metric has changed by an infinitesimal Weyl transformation 

5g^(x) = +2r(x)g^(x). (3.19) 

A similar infinitesimal transformation exists for the more general case ( |3.11| ) 

p -> e~ 2r p ~ (1 - 2r)p, x^ 1 -> + a/*, (3.20) 

a " = \f dfJP(ff)gr{x,f/)d v T(x), (3.21) 
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The effect of this transformation over the leading term of the metric is again (jUSD . The 
leading terms of the scalar field transforms as 



-Auv) = _( d _ & uv ) T (x)^ d - Auv \ 5^ d ~ Al ^ = -(d - A IR )T(x)^ d ~ AlR \ (3.22) 

Subleading terms in the expansion of the metric and the scalar field ( |3.14| ) have a covariant 
dependence on the leading terms, so the effect of this transformation is to make them 
dependent on the combinations gffi + 5gffl(x), (j)( d ~ Auv ^ + 5(j)^ d ~ Auv h There is however 
one term that does not change covariantly, due to the logarithmic term in the expansion 
of the metric ( [3. 14 ). We have then, that when p — > 0, 



Sa (d) _ ^^Sa {0) + 9 ^ 56 (d ~ Auv) - 2r(x)hW 



(3.23) 



while, when p — > oo, 



* = w g « + w^ 6 ^ d ~ AlR) - 2r(x)/ * (3 - 24) 

Introducing a spurion field requires to promote the infinitesimal PBH transformation 
to a finite transformation, including non-linear effects 

p e' 2r p + 5p, x^^x^ + a^, (3.25) 

where 

i rf e ~ 2T 

dp'i 2 {p')g^(x,p')d u T{x) + 5a^. (3.26) 



2 

The boundary values of the fields are affected by this transformation, so in general the 
resulting solution will not interpolate between two AdS spaces any more. However, if the 
derivatives of r(x) are small enough we can apply the same derivative expansion ( [3.14 ) in 
the near-horizon region up to some value of the radial coordinate pm. 

Indeed, the non-linear generalization of the PBH transformation is consistent with an 
expansion of 5p and Sa^ in terms of the number of derivatives of r 



5p = 5p^(d 2 T) + 5p( 4 )(S*T) + 5p (6) (d 6 r) + • • • 
Each term is fixed by demanding that the metric keeps the form 



5a» = Sa» (d 3 r) + Sat, (<9 5 r) + 5a» (d 7 r) + • • • (3.27) 



4p 2 p 



ds 2 = t{pe- 2T )-f^ + -e 2T g llu (pe- 2T ,x)dx fi dx u . (3.28) 



For instance, the contribution to the shift in p that is second order in derivatives is deter- 
mined by the equation 



p ^ I pe~ 2T 



SpW = -EZJLg°0d a TdBT. (3.29) 
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Although straightforward, this expansion can be quite cumbersome. It is also not com- 
pletely unambiguous, for instance the equation above has a homogeneous solution 

6$ =pe- 2T ~f(x), (3.30) 

for some arbitrary function f(x) that has to be second order in derivatives in order to 
be consistent with the expansion, an example is f(x) = e~ 2r g^d a TdpT . Close to the 
boundary, this can be seen as a field redefinition of t adding higher derivative terms 

r -)• r + J—e-^g^daTdpr + 0(8 A ). (3.31) 

The final effect of the transformation is to do a finite Weyl rescaling of the sources for 
the metric and the scalar field, both at the boundary and in the near-horizon region 

0(d-A /fl ) _^ e -{d-A IR )r{x)^d-A IR ) _ £ ir _ ( 3 _ 32 ) 

All other terms in the derivative expansion depend on these sources and transform covari- 
antly according to the Weyl transformation, except for the anomalous contribution 

9® [<$ , (d " A) ] -+ e-^W (<#> [g, u , 4>] - 2r{x)h$ [g, v , 0]) (3.33) 

The additional contribution appears only in even d dimensions, and is due to the logarithmic 
term in the expansion of the metric in ( p. 14 ) . 



3.3 Anomalous contribution to the generating functional 

Having established how the spurion should be introduced in the holographic dual, we can 
now derive the anomalous contribution to the generating functional and check if the results 
agree with the expectation from field theory. 

According to the AdS/CFT correspondence, the generating functional of the field 
theory is the on-shell action of the gravitational theory, as a function of the boundary values 
of the fields. We will assume a large- N, strong coupling approximation, such that the on- 
shell action reduces to the classical gravitational action fl2,3|) evaluated on the solutions. 
We will define the on-shell action with cutoffs in the radial direction, puv "C 1 an d pm 3> 1 

1 r fpm 

S[puv,Pir] = -TTli / ddx / dpV^G (R + 3 M 4>d M 4> + 2V (</>)) 

lK J J orrv 

(3.34) 



+ 4r / <r'.r\'-GK 
K 



PUV 

1 



PIR 



d d xV-GK 



Puv 



On the second line we have added Gibbons-Hawking terms, since we will do variations of 
the classical solution that affect the value of the metric at both cutoffs. Here K is the 
extrinsic curvature, defined from the metric components G^ u = hg^u- \/ —G is the square 
root of the determinant of G^ u , while y— G is the square root of the determinant of the full 
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metric. N 2 = g pp = is the lapse function. The explicit form of the extrinsic curvature 
is 

K„ u = T^dpG^ = j(d p - ~\ g pu , K = G^K^ = pgTK^. (3.35) 



In anti-de Sitter space, on-shell action has the form [32,35] 

L 



S ^ 2k 2 



Puv 



dpd d x p-i- 1 b(x, p) , (3.36) 

where, close to the boundary p ~ puy — > 0, 

oo 

b(x,p) = Y J bn(x)p n + -- - . (3.37) 

n=0 

Where the dots refer to possible contributions that are not integer powers of p, which can 
appear in the presence of matter fields in the bulk. From this expansion one can compute 
the form of the divergent contributions to the on-shell action. Note that in even dimensions 
there is a divergent logarithmic contribution from the order d term, 

S div D-^ J d d x log(p uv )b d (x), (3.38) 

that in fact is proportional to the anomaly 

b d = - d -L d - 2 A d . (3.39) 

Where here and in the following the anomaly polynomial A d has been defined in such a 
way that coefficient of the Euler density E d is one. The divergence is eliminated by adding 
a counterterm S c t = — Sdiv 

Under an infinitesimal PBH transformation the p coordinate is rescaled by a e _2r 
factor, so that a change of variables in the integral appearing in ( 3.36j ) puts this factor in 



the limits of integration puy — > pyye 2r , pm — > pme 2t . Then, the on-shell action is 
shifted by a finite term proportional to the anomaly 

$S div = ^- J d d xT{x)^-gi»)A d . (3.40) 

This can be compensated by a change of the finite part of the action if one combines the 
PBH transformation with a Weyl transformation of the boundary metric, in such a way 
that the boundary metric is left invariant. Note that this term is scheme-dependent, in 
a different scheme one could choose the cutoffs to depend on the spacetime coordinates, 
puy — > puve~ 2a an d pm — > pme~ 2a , in such a way that the value of r is shifted r — > 
t + a. Demanding that there are no anomalous contributions before performing the PBH 
transformation will fix the scheme to a = 0. 
In our case the action takes the form 

S = -4 / dpd d x£(p)p-2- 1 ^- g mb(x,p). (3.41) 

K J PUV 
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We can expand the function b(x, p) in the AdS regions 

oo 

p~puv, b(x,p) = Y,b U n V {x)p n , (3.42) 

n=0 

oo 

p~p IR , b(x,p) = Y,b I n R (x) P n . (3.43) 

n=0 

The UV expansion is the usual p — > expansion. The IR expansion is a derivative expansion 
that is valid for finite p deep in the IR AdS region but not all the way to the horizon, where 
the expansion breaks down. The coefficients b% v (x) and b^ R (x) take the same form as 
( 3.39 ), with L changed to Ljjy and L/_r, respectively. The factor £(p) becomes a constant, 
Ljjy or Lir, plus /3-dependent corrections. Since the flow is driven by the matter fields, 
these contributions take the form of non-integer powers of p. 

From these expansions we can infer that there are contributions proportional to the 
anomaly coming from both the UV and IR cutoff. 

SD ~^I ddx { LlRd ' 1 l ° g{PlR) ~ Luyd ~ l log ( puv) ) Ad ' (3 - 44) 

Following the same arguments as for the pure AdS case, a PBH transformation will add 
to the action a term linear in r and proportional to the difference between the UV and IR 
anomalies 



d 



$S = -^ [L uv d - 1 - Lir j I d d xT{x)^-g^)A d . (3.45) 

The contribution to the energy-momentum tensor of the theory can be computed from a 
variation of the expression above with respect to the metric. We will now give an alternative 
derivation that will serve as a cross-check of our result. Using the equations of motion, the 
on-shell action can also be expressed as 

S[puv,Pm] =-i>J ddx d P^G Qtt^G^ + H d p( ^j . (3.46) 

This has the usual form of a classical action with a vanishing Hamiltonian, a well-known 
result in the ADM formalism that has also been used in the Hamilton-Jacobi formulation 
of RG flows [26,33]. The conjugate momenta to the metric and the scalar field are defined 

as 

"" SS K^-KCT, ^= K *^ = dp <t>. (3.47) 



SS[p uv ,pm] = -^2 J d^xy/^f^SG^+^sA 



Because the action is on-shell, a variation of the action only has contributions from the 
cutoffs 

Puv 

(3.48) 

PIR 

The variation at the boundary should vanish for dynamical modes, this corresponds to the 
Dirichlet boundary condition 5(f) = 0, 5G^ V = 0. 

The regulated expectation value of the energy-momentum tensor is 

= --V^TV- (3-49) 
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At the boundary puv 0, this is divergent and one needs to add counterterms that cancel 
those divergences, which can be done by adding covariant local terms at the cutoff [29]. In 
the near-horizon region the regulated energy-momentum tensor is finite and no additional 
contributions should be added. 

After eliminating the divergences, the puy — > limit leaves only finite contributions to 
the energy-momentum tensor. After the holographic renormalization procedure has been 
implemented, the variation of the action takes the form 



5S = J d d x^gV) [~ (V) 5g$ + (O) 8^™) 



(3.50) 



In our case the expectation value of the energy-momentum tensor has contributions from 
the boundary and the near-horizon region. Near-horizon contributions depend in general 
on pm, that acts as an IR regulator. The exception is the order d term (with no logarithmic 
factors), that has the same form as the renormalized boundary contribution. Using the 
results of [29,35], the contributions of order d are 

W d b (0) , 4> {d ~ Auv) ] = ^ ( W 1 " 1 - W- 1 ) $ fo<°> , ^-)] + • • • (3.51) 

where the additional terms depend on lower order terms g( n<d ) or gffjg^ a ^9^p- 

When we do the non-linear PBH transformation, the energy-momentum tensor be- 
comes dependent on the dressed sources ( 3.32j ). The transformation is fully covariant, 



except for the additional shift of the order d component ( |3.33[) , that produces an addi- 
tional contribution to the energy-momentum tensor 

S(T^) [gj uv ] = -t(x)± (W" 1 - L IR d - 1 ) e^-^h^g^hv]- (3.52) 

As was explained in [29] , tiff equals the variation with respect to the metric of the confor- 
mal anomaly, so this implies that the generating functional contains a new term 

r d ^ (L uv d - 1 - L IR d - 1 ) j d d x V /Z ^r(x) A d , (3.53) 



In agreement with our previous derivation ( 3.45| ). Note however that the variation of the 



Euler density is zero, so strictly speaking this derivation only captures the contributions 
to the conformal anomaly that are Weyl- invariant. 

The overall coefficient of the Euler density in the new term is proportional to the 
difference of the a-central charges of the UV and IR fixed points [35] 

a uv -a IR x^ (L uv d - 1 - Llr*- 1 ) . (3.54) 

The holographic a-theorem of [10] implies that Ljjv ^ Lir, so indeed the coefficient of 
the anomalous term is positive. For dimensions larger than four this is a prediction from 
holography, there is no proof in field theory of this statement. The holographic calculation 
may be a hint that such a proof actually exists. 
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For two or four dimensions the result is in agreement with the field theory analysis, 
but note that the trace of the new contribution to the energy-momentum tensor vanishes 
g (0) iwg (T^) = 0, so the new contribution does not affect to the trace anomaly in the one- 
point function. In spite of this, in our calculation the trace anomaly is also proportional 
to the difference of central charges, but it is determined by the trace of the g$ term 

<T (T^) [9, 4>uv] = 4r {L uv d ~ 1 - Ljr*- 1 ) A d . (3.55) 

This is the expected result based on physical grounds, by introducing the cutoff pm we are 
discarding the massless degrees of freedom at the IR fixed point, so the anomaly generated 
by degrees of freedom above the cutoff is proportional to the total anomaly minus the 
anomaly of the IR CFT. 

4. Spontaneous breaking and Goldstone boson 

When conformal symmetry is spontaneously broken we expect in principle that a massless 
Goldstone boson associated to dilatations will be present at low energies. 5 According to 
the anomaly matching arguments of [1] its effective action should contain an anomalous 
term similar to the one we have computed for the spurion field. Although it is encouranging 
that one should find the same result for the generating functional as in the field theory 
calculation, there are some caveats in the use of the spurion field that make a separate 
analysis of the dilaton mode necessary. The caveats are the use of a radial cutoff in order 
to compute the generating functional in a derivative expansion and the fact that introducing 
the spurion field alters the boundary values of the metric and other fields; according to the 
usual dictionary of the AdS/CFT correspondence, the spurion should not be interpreted 
as a dynamical mode of the dual field theory. 

Our goal in this section is to study correlation functions of the energy-momentum 
tensor and the scalar operator whose expectation value breaks conformal invariance. When 
the symmetry is broken spontaneously the behaviour of these correlators at low momenta 
is constrained. In particular, a massless pole should appear in the transverse part of mixed 
tensor-scalar correlator, due to the existence of zero energy states. 

In order to compute correlation functions we need to study small fluctuations of the 
metric and the scalar field. Related works are the calculation of correlation functions 
[15, 25-27] in a geometry dual of a particular state in the Coulomb branch of N = 4 
SYM [8], that is an example where there is a massless pole that one can identify with 
the dilaton. In this case the near-horizon geometry is not AdS, but a singular geometry 
and there is a mass gap in the spectrum. In [36] a calculation in a four-dimensional flow 
between two fixed points was carried out, with the result that at small momenta the two- 
point function of a probe scalar becomes that of the IR CFT. A calculation of correlators 
in a two dimensional flow was also made in [37] . 

Although for a general superpotential we were not able to find the exact solutions to 
the equations of motion, we can do an expansion in momentum that is possible to solve 

5 In d — 2 dimensions the analysis becomes more subtle, since IR divergences prevent the appearance of 
free massless modes, although it is still possible to have massless scalar states in the theory [34] . 
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order by order. In order to capture the low momentum behavior of the correlation functions 
this expansion is sufficient. The expansion however breaks down very close to the horizon, 
where regularity or other physical conditions have to be imposed in order to fix the solution. 
Fortunately, we can find the near-horizon solution using a different approximation, that in 
this case is to approximate the near-horizon geometry by AdS. There is a region in the 
geometry where the two approximations overlap, and the matching between the two fixes 
the form of the correlators at small momenta. In appendix we check that this procedure 
indeed captures the leading low momentum behavior in a geometry dual of N = 4 SYM in 
the Coulomb branch, where the full solution can be computed analytically. 

4.1 Metric and scalar fluctuations 

In the models we are considering the dual of the RG flow is described by Einstein gravity 
coupled to a scalar field. The equations of motion are 

Rmn = -d M 4>dN<t> - -j^-r 9 mnV (((>), V 2 4> = (4.1) 
a — 1 o<p 

We will adopt Gaussian coordinates for the metric, 

ds 2 = g MN dx M dx N = dr 2 + g^(r, x)dx fi dx v . (4.2) 

The choice g^ r = 0, g rr = 1 corresponds to fixing d+1 diffeomorphisms up to transfomations 
of the form 

<W = V (Af 6v) = 0, if M, N = r. (4.3) 

For the background metric g^ v = e^^'r}^, we can distinguish between d-dimensional 
diffeomorphisms 

£, = 0, e„ = e 2A ^a,(x), (4.4) 
and translations in the radial direction 

Cr = -ff(x), £m = d^(x)e 2A ^ ^ dr'e- 2A ^'\ (4.5) 

For a a constant, these transformations become scale transformation in the AdS back- 
grounds, so it is natural to associate them with dilatations. 

We will expand the equations to linear order in fluctuations 

9»v = e 2 (r/^ + hfj, v ), (f) = (f) + ip, (4.6) 

and use the flat metric to raise and contract indices. For instance, tha trace of the metric 
fluctuation is defined as h = h% = rj^h^. We will use an expansion in Fourier modes. 
Each mode of the metric admits the following decomposition in transverse and longitudinal 
parts, traceless and not traceless 

v = C T + C L + C + ^- ( 4 - 7 ) 
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More explicitly, we define the projectors 



U ^*f3 = p^puP _ _}_p^p*^ ( 4 _ 



so that 



uTT _ 1 tt a/3i 



^ T — 1 P- P^k 

V ~ (rf _ 1)^2)2^ t* v T 

For convenience, we will also define h = h>T + h^, where 

/ lT = ^ = ^^. h L = rrh%, = ^Pfh aP . (4.10) 

The details of the derivation of the equations of motion can be found in the Appendix [A|. 
Here we quote the results for the two modes we will analyze in the following: 

• Tensor fluctuation: 

o=hir+dA'hir-Q 2 e- 2 xi 

• Scalar fluctuations: 

There are two constraints for the trace components of the metric 

= (d - l)A'ti + 2dVcp - 2<))' ip' - q 2 e- 2A h T , (4.12) 
= h' T + 2<j)' (p, (4.13) 

and, defining h! = e~ 2A H, two second order equations for the metric and the scalar 
fluctuations: 

= e~ 2A H' + j^jdVp + WW, ( 4 - 14 ) 
= <p" + dA'tp' - d 2 Vip - q 2 e~ 2 \ + ^-e- 2A H. (4.15) 

4.2 Tensor fluctuation 



TT 



We start with the simpler case of a tensor fluctuation, given by ( 4.11 ), with h 
£^ u (q)H T . We will also define Q = \j—q 2 = \/uj 2 — q 2 . 

In general, we cannot solve ( 4.11 ) exactly but we can solve it perturbatively for small 

Q 2 

= + Hf J dre' dA + 0(Q 2 ), (4.16) 
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where Hq and Hf do not depend on the radial coordinate but may depend on q. One can 
solve recursively for the next orders 



H T = Y J Q 2n H {n \ 



(4.17) 



n>0 



where H^ = and the recursive equation is 

hW" + dA'H^' = - e - 2A Q 2 H {n -^ 
The solution to this equation takes the form 

H (n) = _ j dre -dA j dre (d-2)A Q 2 H (n-l)_ 



(4.18) 



(4.19) 



One should actually keep terms up to order 0{Q d ) (for d even), since the normalizable 
term is of this order. 

If Q 2 is small enough this solution can be a good approximation even in the near- 
horizon region, assuming Q 2 e~ 2r l LlR <C 1, 



1 + 0{Q 2 e- 2r ' LlR ) 



+ e 



-dr/L IR 



LlR 

~d~ 



(4.20) 



This expression is valid for d > 2 and even. The perturbative solution breaks down very 
close to the horizon at r — > -co since the factor e~ 2A that multiplies Q 2 in the equation 
becomes dominant. Therefore we cannot use this solution in order to impose boundary 
conditions in the IR. For that we will have to solve the equation of motion in the near 
horizon region, where it takes the form 



= H T " + 



d 
Lir 



H T> _ e -2r/L IR 2 H T_ 



(4.21) 



The solutions to this equation are Bessel functions. There are several possible choices of 
boundary conditions at the horizon. We impose ingoing boundary conditions at the hori- 
zon, following the prescription of [38]. This condition basically means that the solution 
behaves as an ingoing plane-wave in the vicinity of the horizon. Note that if the solu- 
tion is continued analytically to Euclidean signature, this choice corresponds to fixing the 
exponential behavior of the solution, in such a way that the solution is regular. 

The solution satisfying ingoing boundary conditions at the horizon takes the form 



Hl = C [L IR Qe- r / L ^) d/2 Hf (L IR Qe^ 



•/Li 



(4.22) 



where H^ is the Hankel function of the first kind. 

In the low momentum limit we expect that there is a region in the geometry where 
both the solution HT and the solution HT are valid, this will be when r — > — oo but 
L IR Qe- r / LlR < 1. In this case we should be able to match the leading contributions of 
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both solutions ( 4.20| ) and ( 4.22 ). In the case at hand, for d even dimensions, the solution 
(4.20) in the limit r — > — oo takes the form 



H T _^ H T e -dr/L IR + ... 



while the solution (|42|) in the limit L m Qe- r l LlR «; 1 is 



Hi ~C 



+ • • • + ^(^Q) d log(L /H Q) 2 e-^ i - + 

7T 7T 



(4.23) 



(4.24) 



where for d = 2,4,6, = 2,4,6 and bd = 1/2,1/8,1/48. There are analytic terms 
proportional to (Lj R Q) d , as the term proportional to Hq in (|4.20 ), but for d > 4 those 
will only introduce contact terms in the correlation function and we can ignore them. We 
will comment on the d = 2 case when we compute the correlation function. 

In both cases we have identified the leading terms of the normalizable and non- 
normalizable solutions. Note that the non-normalizable solution can have contributions 
that are less suppressed than the normalizable solution when r — > oo. The matching of the 
two solutions determines 

d b d 
n n > n-i — ~ 

ad 



-{L IR Q) d \og{L IR Q) 2 H { 



Lir ad 

Then, the ration between the non-normalizable and normalizable solutions is 

H? d b d . 



(4.25) 



Hi 



-{L IR Qf\og{L IR QY 



(4.26) 



Li R ad 

The dependence on the momentum is in agreement with the two-point function of the 
energy-momentum tensor in a CFT. 

In terms of the boundary value of the metric hb^, we can write the fluctuation as 

1 



h 



TT 



-U..,^h baP 



1 + G TT {q)e -dr/L uv + 



where, at low momentum 

G TT { q ) 

4.3 Scalar fluctuation 



Lyy b d 
Li R ad 



(L IR Q) d log(L IR Qf 



(4.27) 



(4.28) 



A scalar fluctuation involves both the scalar field and the scalar components of the metric, 
that are coupled through ( (4.14 , 4.15 ). 

We can now decouple the scalar field by multiplying ( [4.15 ) by e 2A /(f>' , taking the 
derivative once with respect to the radial coordinate, substituting H' by the expression 
obtained from ( [4.14 ) and multiplying by 4>b e ~ 2A - The result is a third order equation for 
(p, for which there are three solutions. The form of the equation is quite involved, but we 
can simplify things by factorizing out one solution. This is possible thanks to the residual 
diffeomorphism (45), under which the scalar field transforms as 

5(f) = (f)' £ r = -<r{x)dW. (4.29) 
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where d is a derivative with respect to the background field cpQ. The metric is also affected 
by this transformation 



5g»u = d^ u + dviii + Cr-drSV = 2e 2A (d^d v a{x) j dr't 



-2A 



W 



d 



-Vnuv(x)). (4.30) 



The gauge mode ip = dW is a solution to the equation of motion. This fact enables us 
to factorize it as follows 



where 



a 



(d 2 + ai d + a )(d 



t d 2 W (d + 2)W 



d 2 W \ 
dW J 



<p = 0, 



(4.31) 



dW (d-l)dW 
Ad 2 WW 



q 2 e 2A 



+ 



2dW 2 



dW 2 (d - l)dW 2 (d - lfdW 2 



2d 3 W 
dW 



2. 



(4.32) 



Note that the gauge mode is a valid solution for any value of the momentum and 
of the radial coordinate and that it vanishes at the boundary in the way expected for a 
normalizable solution ip ~ <pQ. The transformation of the metric ( 4.3CQ on the other hand 
does not vanish at the boundary. Then, imposing a Dirichlet boundary condition on the 
metric will fix this mode. 

We proceed to find the two other solutions. In order to simplify things further, we 
define 



d 



d 2 W \ 
dW ) 



W 



(dwy 

Then, G satisfies the following second order equation 

-2A 



: e- dA G. 



d z G + dBdG - <f 



{dwy 



rG = 0, 



(4.33) 



(4.34) 



where 



B = 2 \og{W) - \og{dW) - (d - 2) A. (4.35) 

In general, we cannot solve ( [4.34 ) exactly but we can use the matching procedure in the 
same way as for the tensor components of the metric. 

First we solve (4.34) perturbatively for small q 2 . The solution takes the form 



G = C x 
+ Co 



1 + 



9oe 



-B 



,B-2A 



#0 



(dW) 2 



boe~ B + q 2 



-B 



+ 0(q 4 ) 

e B - 2A 
(dW) 2 



b e~ B + 0(q 4 



(4.36) 



Then we solve ( [4.34 ) in the near horizon region, where it takes the form 



d 2 G + ^-^dG + ( ±10) 2 - ^) 2 ~- 2 G = 



V A J 



(4.37) 
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The solutions to this equation are Bessel functions. We will pick up a solution satisfying 
ingoing boundary conditions at the horizon (regularity in the Euclidean) 



G = C h [L IR Q((f>, 



H 



(i) 



where 



\ = d — A 



IR 



(4.38) 



(4.39) 



Note that the relation between ip and G ( 4.33 ) is such that close to the horizon <p has 
different power-like dependence on (</>o — (f) m ), but the exponential behavior is the same. 
Then, the ingoing condition on G is equivalent to an ingoing condition for <p. Note that the 
gauge mode does not contribute to the exponential behavior since dW ~ (tfio — (j) m ), so it is 
not affected by the ingoing boundary condition. One can systematically find corrections to 
this solution by expanding the coefficients of the differential equation in powers of (0 — <fi m ) 
and solving order by order. 

Expanding the perturbative solution ( |4.36| ) around the horizon 0o = 4>m we have 



G ~ Ci + C 



L m \ 2 



'(d- l)2(2A-d + 2) ^ u ^' 
Expanding the near horizon solution ( [4.38| ) in small momenta 



2\-d+2 
\ A 



(4.40) 



G~C h 



a + b [L IR Q((f> - 4> m )~ 



2X-d+2 



where 



b 
a 



2 2u T(u)T(u + 1) 



(4.41) 



(4.42) 



7r(l + COt(7TZ/)) 

As we show in Appendix [b], higher order corrections to the near-horizon solution do not 
affect to the leading terms in the low momentum expansion. Matching the two solutions 
fixes 

Go b ^2\-d+2 b ^d~2A IR +2 

-{LirQ) =- 11 



(LirQT 



(4.43) 



Near the boundary the solution for G takes the same form as fl4.40p with the replace- 
ments Lj R — > Luv, A = Ajjv an d 0o — 4>m — > 4>o- Solving ( |4.33| ) for tp and including the 
gauge solution p = C g dW, we can then write its asymptotic form 



ip = fa 



C n ^ + 



Co 



A 



uv 



9 L 



uv 



CiL 



(L UV Q) 2 2(d - l)(2A uv -d + 2) 
d-1 



LuvQ4>o 



uv- 



A uv (d-2A UV ) 

In the following we will define the coefficients of the leading order expansion as: 



(4.44) 



<P ~ Cg^O + C N ^ + C NN ^ 



(4.45) 
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The metric fluctuations h and hr can be computed from ( 4.14 ) and ( [4,13| ) respectively. 
The expansion of the transverse component close to the boundary is, to leading order 



h T = h Tb - 2C g (W - W(0)) - 2 



~ h Tb - C, 



Cn 



2+- 



i (v> - c g dw) 

2A UV 



d 



-Cnn^o 



(4.46) 



A uv + r» 

Here hxb 1S the boundary value of the metric, that we fix by imposing Dirichlet boundary 
conditions. This determines the lower limit of integration, that should be such that no 
additional contributions change the value of h^ b . 
The expansion of the full trace is 



h~h b - - — T C g (pQ - C N - 



+ 1) 



2+- 



d-i-»™ (d-i)(Aw + iy° d{d-i) 

Plugging (|446| ) and (|4i47| ) into ( ^T2|) , we find the following condition 



- + AAuvJAuv -d) ^^ (447) 



C g = {d- 2A 



uv 



(QLuv 



\2' 



(4.48) 



This relation implies that the close to the boundary the solution for the scalar field has the 
same radial dependence as the solutions in the decoupled system, where the background 
scalar is set to zero and the geometry is pure AdS. 



In terms of C g : 



Luv 



Co 



T 2 



(d-1) {L uv Qf 



(d — l)Luv a 



(LirQ) 



d-2A IR 



C-i 



(4.49) 



If we define the coefficients of the boundary expansion of the scalar field and the scalar 
component of the metric as 



we find that 



and 



h = h b + h N <f>* uv ■ ■ ■ 
L uv (d-1) 



where 



G s (q) 



¥b = Ci 

¥N = Cg 

b 1 



A uv (d-2A UV ) 
Auv 



L 



uv 



G s (q)ip b 



Lir 



a 2Ajjy — d Ljjv 
For the metric we find, using ( 4.13 ) and ( A.23| ) 

h NT = G T (q)ip b 
h-NL = G L (q)tp b 



(LirQ) 



d-2A IR 



(4.50) 
(4.51) 

(4.52) 

(4.53) 

(4.54) 



(4.55) 
(4.56) 
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where 



G T {q) 
G L (q) 



d ■ 

2A[/y(2A 



uv 



d-l) 



d(d-l) 



(4.57) 
(4.58) 



4.4 One and two-point functions 



In order to compute the correlation functions of the scalar operator and the energy- 
momentum tensor, we need to evalue the on-shell action for the linearized fluctuations. 
We can use the results of [25,26,31] adapted to our case. The on-shell action to quadratic 
order in the fluctuations is 



(4.59) 



where the variation of the on-shell action coincides with ( p. 46 ) but now the only contri- 
bution is evaluated at the boundary and we are using the Gaussian coordinate r, so the 
extrinsic curvature is = \d r g^ u and the variation is 



5S[g^ v ,4>o\e 2A h^ v ,ip\ = lim — / d d x^g 

r— >co k, z J 



-(K^-Kg^)e 2A h^ + ^ 



(4.60) 



We will implicitly assume the r — > oo limit. In order for the action to be well defined we 
will have to add counterterms that cancel the divergent pieces. The counterterms take the 
form [26] 6 



1 



d d x. 



-9 



W + 



C2 



+ 



(4.61) 



For our analysis it is enough to consider counterterms involving no derivatives. In order 
to compute corrections that are of higher order in the low momentum expansion, we will 
need to add more counterterms depending on the curvature of the background metric and 
derivatives of the scalar field. 

To linear order in the fluctuations, the finite contributions in the r 
lently, the <j)$ — ¥ limit are 



oo or, equiva- 



S(i) 



C2 



d d xe dA 



C-2 



d d x<fb. 



(4.62) 



This shows that the one-point function of the scalar field is a nonzero constant. Here ipb is 
dimensionless (c2 should have dimension one) , while the product of the source for the scalar 
field Jb times the expectation value (O) should have dimension d. Then, ipb is related to the 



6 This choice of counterterms corresponds to a particular renormalization scheme, as explained in [26]. In 
this case the energy density of the vacuum will be zero for the particular flows we are considering. In general 
any function which is a solution to the differential equation that relates the potential to the superpotential 



in (2.4) can be used as counterterm. We want to thank Ioannis Papadimitriou for explaining this point to 
us. 
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actual source by a factor ipb = M Auv d Jb, where M is some mass scale. The expectation 
value of the scalar operator is 



(O) 



6SW 



(4.63) 



We can use this expression to fix the relation between the source Jb and the boundary 
value of the scalar field (fb as 



--(0)J b . 

C2 



(4.64) 



The on-shell action to quadratic order is, in our gauge 

S (2) =-jsj ddx ^9 \h^'h pv - h^'K + ~<p'<p + \{C2 + d 2 W) <p 2 



(4.65) 



where indices are raised with the flat metric rju V . Using the background scalar field (fro as 
radial coordinate, the on-shell action becomes 



S {2) = -^2 J d d xe dA dW 
Expanding close to the boundary as 

—d 1 

f = ^b(t>Q uv + <PN<l>a H , hfu, = h b ^ + h N ^ 

The divergent contributions to the action are 
1 



(4.66) 



(4.67) 



S, 



(2) 



o 



o 00 <Pb + o ^0 Vb 



This imposes the condition 
Then, (|4.64| ) becomes 



C2Luv = d — 2A 



uv- 



k 2 L 



P>b 



uv 



2A.UV - d 

The finite contributions to the action are 



(O) J b . 



(4.68) 



(4.69) 



S, 



1 



(2) 



K 2 L 



uv 



d d x 



d 



d, v d-2A uv 2d-3A uv 

N ^ h bu + j ^ NLfb + 4 h ^ b 

(4.70) 



htf hb^u + sh N u hbv + 
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Using the decomposition fl4.7| ) and ( [4,9|) and assuming hjf 1 



0, we can write the 



action as 



S, 



(2) 



d d q 

k 2 L uv J (>p 



TT 



1 



b a/3 



+h 



Nfw („2\2 



(q 2 ) 2 (d 

d a d-2A uv 2d-3A uv a 

+ g h N ah u + g VNVb + | K 

The solution we have found depend on the boundary sources as 

1 



■TT 



h 

h NT 



(q 2 ) 2 
1 



n^G TT ( q )h b a 



/3' 



P^hNfiv = G T (q)(f h , 



1 



hNL = -^Pl hNnv = G (q)<pb, 
ifN = G s (q)tp b . 



1 TD a P 



Therefore, using that \P T + -^P L 



5, 



(2) 



k 2 Ljjv J (2vr) 



d d q 



1 



'^(q 2 ) 2 
^PfG T s {q) + ±PfG L s {q) 



m v ^G ll {q)h b a0 
d - 2A UV 

lib aj3 H 7, 



ip b G s (q)(p b 



where 



Therefore we find 



G T s [q) = d -^ [Z ^G T (q) + dG L (q) + 2(2d - 3A UV ) 
G L s {q) = dG T {q) + 2{2d - 3A UV ). 



G T (q) = 2A uv (2A U v-d) +G ^ 
G^q)=4(d-2A UV ), 



G T s {q)-G L s {q) 



2A C /v(2A C /y - d) 
d-1 



(4.71) 

(4.72) 

(4.73) 

(4.74) 
(4.75) 



(4.76) 

(4.77) 
(4.78) 

(4.79) 
(4.80) 
(4.81) 



The two-point function of the energy-momentum tensor can be obtained by taking the 
derivative of the action twice with respect to the boundary metric, and at low momentum 
is determined by (|1.28| ) 

L IR d 1 



T^r a n ~ h d 



k 2 L UV (q 2 ) 2 



U^Q d log(L IR Qy. 



(4.82) 



where hd is a constant depending on the number of dimensions. Note also that the overall 
coefficient is dimensionless. This is the expected behavior when the dual field theory flows 
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to an IR CFT. As expected for spontaneous breaking of conformal symmetry, there are no 
contributions to the trace of the energy-momentum tensor in the two-point function. 

The scalar correlation function is obtained from ( 4.76| ) taking the derivative twice with 
respect to the source J b , defined in ( p~69l) , 



K 2 Tr^r K 2 T r d+2-2A IR 

(OO) = " LUV . (O) 2 G s (q) ~ s d ^- {Of Q rf " 2A -. (4.83) 

Here is a numerical constant depending on the number of dimensions d and the dimension 
of the scalar operator Ajjy. Note that the dimension of the two-point function agrees with 
its UV value 2Ajjy — d. Since Am > d, there is a singularity at small momentum, which 
is in agreement with the expectations from field theory. However, if the dilaton was a free 
massless mode, we would expect this singularity to be a pole, we will comment more on 
this in the Discussion. Note that for a generic value of Am a pole cannot appear at higher 
order in the momentum expansion because the overall singular factor will be multiplied by 
integer powers of Q 2 . On the other hand, if Am is an integer such a pole could appear, 
but checking if this is the case would require a separate analysis, since the expansion we 
have used throughout the paper assume that Am is not an integer. 

The mixed correlator of the energy momentum-tensor and the scalar operator is ob- 
tained by taking a derivative with respect to each source: 



Lp^ G T {q) + l_p^ G L {q) 



9 T 8 VlJ 1 9 



4{2A UV - d) 

Writing the projection operators explicitly we are left with 

^P^(G T s -G L s ) + t v G L s {q) 



(4.84) 



q 2 ' 



4{2A UV - d) 

-W^7'*'-< >'>~. (4.85) 

The first term is the one expected when conformal symmetry is spontaneously broken. The 
second term is a contact term, that appears in the retarded correlator 

(T»»(x)0(y)) R ~ (0)rrsW(x-y), (4.86) 

note that it is a real term, so it does not appear in the commutator. A similar term appears 
in the correlation function computed in the Coulomb branch flow [25,26,31]. In [31] it was 
argued that this term appears because the generating functional has a term of the form 
~ (O)j Jb, which gives an additional contribution. Substracting it one would recover the 
true energy-momentum tensor of the field theory in the absence of sources. 

We have argued that the spurion does not correspond to the Goldstone boson of spon- 
taneous breaking of conformal invariance according to the usual holographic dictionary 
between normalizable solutions and dynamical modes. However, correlators of the energy- 
momentum tensor and scalar operator have a low energy-momentum behavior which is 
in agreement with the expectations from field theory. So what is the dilaton? It should 
correspond to a normalizable solution of the equations of motion at zero momentum. The 
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normalizable solution in ( 4.36| ) is the one with C\ = 0, but the condition ( 4.48 ) at q 2 = 
implies that Co = as well. Therefore, the only possible solution has the same form as the 
linearized form of the spurion, but with appropriate boundary conditions for the metric: 



ip = T (x)dW, V = -2t(x)(W -W(0)) VfMU . (4.87) 

Note the similarity with the linearized form of the spurion, the difference is in the boundary 
value of the metric, that for the dilaton mode is fixed while for the spurion it changes by 
a Weyl factor. 



4.5 Comparison with field theory 

Let us first show that the results we have obtained are consistent with low-energy theorems 
for a field theory with spontaneous breaking of symmetry. There are two versions of the 
theorem, that can be found in textbooks, e.g. [39] and that we will skecth again here for 
convenience. The first version is based on the effective action. Suppose the action and 
measure of the path integral are invariant under a continous symmetry, and for simplicity 
let us assume translation invariance and constant fied configurations 

5(f) m = iet mn (f) n , (4.88) 

where (f) n are generic fields, e a small parameter and t mn the symmetry generators. If the 
transformation is a dilatation then t mn is the matrix of conformal dimensions. The effective 
potential will be invariant under this symmetry 

TX^ Wn = 0. (4.89) 
Taking a derivative of this expression with respect to <pi , one finds the following condition 

In a vacuum configuration, the first term vanishes since it corresponds to a minimum of 
the effective potential. The second term can be identified with the inverse of the zero- 
momentum propagator 

^^(0)^ = 0. (4.91) 

m,n 

Therefore, there is an eigenvector of the zero-momentum inverse propagator with zero 
eigenvalue. This is usually interpreted as a pole at zero momentum for the propagator 
itself. However, from this argument one can only establish the existence of a singularity 
at zero momentum, not the nature of the singularity itself. The result we find for Q4.83 ) is 
then consistent with this version of Goldstone's theorem. 

The second version of the theorem is based on the Ward identities of the theory. For 
concretness we will specialize to space-time transformations. Using the energy- momentum 
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tensor we can construct the following charges that generate translations, Lorentz transfor- 
mations and dilatations: 

P* = f «- = /AJ^, <«■«> 

When conformal symmetry is spontaneously broken by a constant expectation value of an 
operator O of dimension A, the commutater of the generators with the operator vanish for 
translations and Lorentz transformations, but not for dilatations. 

{[P»,O]) = 0, {[M' w ,O]) = 0, ([D ti ,0)) = -A(0). (4.93) 

On the other hand, since the symmetry is only spontaneously broken, the currents associ- 
ated to the symmetries should be conserved 

dy([T^(y),O(x)}) = 0, dy([y^T^(y),O(x)])=0, W {[y v T v »{y), 0(x)\) = 0. (4.94) 



The next step is to write the commutator between the energy-momentum tensor and the 
scalar field in terms of a spectral decomposition: 

([T^(y),0(x)})= j^yr p^(p)e ip( - y -^ + p f * u (p)e- ip{y - x) . (4.95) 



, (2vr)3 . 

From ( 4,93| ) and ( |4,94j ), one finds that the spectral function should take the form 

(Tip) = (^V - pV)l?(po)/>( V), (4.96) 

with the condition 

p 2 p(p 2 ) = 0, J dp 2 p{p 2 ) = A (O) . (4.97) 

The two conditions can be satisifed only if 

p(p 2 ) = A(0)5(p 2 ). (4.98) 

A delta function in the commutator between T^ u and O implies that there should be a pole 
in the real part of the time-ordered correlators, proportional to A(0). The holographic 
calculation ( fi.85 ) indeed shows it. 



4.6 Promoting the spurion to a dynamical mode 

Since the spurion behaves as a regular and normalizable mode for the scalar field, we 
can promote it to a dynamical mode by introducing a new set of boundary conditions 
for the metric. This may be interesting in the context of effective low energy theories for 
applications to condensed matter, as it can be an example of "holographic deconstruction" 
following the proposal made in [40]. In this case, the interesting IR dynamics is captured 
by the geometry below some cutoff in the radial direction. The IR cutoff effectively gaps 
the modes living between the cutoff and the boundary, but the information about the UV 
physics enters through the coupling of the modes below the cutoff with a massless mode, 
a "Goldstone boson" with a radial profile between the IR cutoff and the boundary. The 
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spurion field will be one of these modes for theories with spontaneous breaking of conformal 
invar iance. 



Starting with the on-shell action (3.34) we can impose Dirichlet conditions for the 
scalar field (we will fix the source to vanish, so this is allowed) and a mixed condition 
for the metric if we add additional boundary terms. First, we split the variation into the 
traceless and the trace contributions 



5G 



0. 



We will impose a Dirichlet condition as usual for the traceless part 5G T 



/it/ 



(4.99) 
0, but for the 



trace part we will introduce a Neumann condition. In order to have a good variational 
principle we should add boundary terms to the action that cancel the contribution from 
the on-shell action S t ot = S + Ijjy + Iir 



d-1 
2n 2 p uv 



d d xJ-G(p 



uvj 



Iir 



d-1 



2k 2 p 



IR 



d d x 



-G(p 



IR) 



(4.100) 



Note that \J—G in ( |3.48| ) has an additional factor £/2p. 
Now, using Q3.48P , the total variation of the action is 



&Stot [puv, Pir] = 5S[p uv ,piR\ + SIuv + Sim = 0. 
Imposing the Dirichlet conditions 5(f) = 0, SG T UV = 0, this takes the form 

1 



$S t ot[puv,PiR] 



2k 2 



d a x 



-G ( G^ v 



d(d-l) 



5t(x) 



Puv 



Pir 



The variation vanish for any 5t{x) if 



g aP pd p g a p 



Puv 



PIR 



(4.101) 



(4.102) 



(4.103) 



Since the variation is evaluated at the cutoffs, we can use the expansion ( 3.14 ) in the bound- 
ary condition ( |4. 1Q3| ) . At the boundary there are divergent terms that should be removed 
with local counterterms. The reminder in the puy — > limit is a finite contribution. For 
d > 2, the leading contribution in the derivative expansion of ( f4.103j ) is 



PIR 9 



(P)ap (2) 



+ ... 



o. 



(4.104) 



This will determine the equations of motion of the dilaton. For a flat background metric 
9av = Vpu, the infinitesimal PBH transformation (3.20) changes the metric to 

g$ ~ Ljj^d^Sr, (4.105) 

so, to leading order, the boundary condition ( 4.103[ ) is satisfied for a massless mode 

D8t ~ 0. (4.106) 



7 Although the Weyl transformation also changes the traceless components, we are ultimately interested in 
the effective action of the dilaton in flat space, where these components vanish and the Weyl transformation 
does not affect them. 
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The full equations of motion can be derived from the effective action of the dilaton, that can 
be constructed by computing the generating functional as before including the boundary 
contributions (4.100) and performing the non- linear PBH transformation. The anomalous 
term with a coefficient proportional to the difference of central charges will also appear in 
this case. 



5. Discussion 



We have shown how the conformal anomaly matching arguments presented in [1, 2] are 
realized in theories with holographic duals. The dilaton/spurion is introduced through a 
coordinate transformation and our analysis shows that the effective action contains a term 
whose variation is the difference between the anomalies of the UV and IR theory. 

8 a T oc a (A u d v - A 1 /) . 

The term is naturally present in even dimensions and absent in odd dimensions, where 
there is no anomaly. 

The analysis we have presented is valid for theories whose holographic dual consists 
of Einstein gravity plus matter fields, with a geometry that interpolates between two AdS 
spaces of different radii. The null energy condition guarantees auv > am and makes the 
coefficient of the anomaloues term positive definite. This is a general statement within the 
present, best understood, class of holographic models, however it is not quite close to be 
a general statement about field theories (even leaving large- N and strong coupling aside). 
Note, in the first place, that in the holographic theories we are considering the coefficients 
of the anomaly polynomial at the fixed points are not independent, but there is a single 
independent coefficient which can be taken to be the coefficient of the Euler density. For 
instance, in four dimensions 

A 4 = cW 2 - aE A , (5.1) 

c = a is fixed. In particular this means that the holographic c-theorem (or rather a- 
theorem) does not apply only to the a coefficient but to all the other coefficients of the 
anomaly polynomial. We know this to be not the case in general (e.g. [41,42]). It will 
then be interesting to extend the analysis to cases where a ^ c. In holographic models 
this can be achieved by including higher derivative curvature corrections, an example is 
Gauss-Bonet gravity in five dimensions, for which a holographic a-theorem also exists for 
the a coefficient [46]. 

Let us now comment on the behavior of the scalar correlator for the case of spontaneous 
breaking, whose low-momentum behavior we have found to be 

(OO) ~ (\/^2) a! - 2A " i . 

where Am > d. Note that this is more singular than the pole behavior we expected 
to see. We can compare our results with the results for the singular Coulomb branch 
geometry. There, it was found that there is a pole in this correlator [25,26,31]. Although 
we understand this as coming from the properties of the solution close to the horizon, we 
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are not sure about the physical origin of this difference in the two cases. A qualitative 
difference is that in the Coulomb branch the spectrum is gapped, which probably limits 
the singular behaviour of the scalar correlator to a pole, while in the cases we study there 
is always a large number of massless degrees of freedom. It would be interesting to explore 
other models with spontaneous breaking of global symmetries, and compare the two-point 
function of the operator that triggers the breaking in cases where the spectrum is gapped 
with those where there are more massless states beside the Goldstone bosons. 

There are several examples where massless degrees of freedom affect the qualitative 
behavior of correlators. For instance, in the Landau-Ginzburg model long-range interac- 
tions can in general modify the momentum dependence of the correlation function of the 
order parameter, leading to non-analytic dependence on momentum [43]. Long-range in- 
teractions can also modify the dispersion relation of Goldstone bosons in non-relativistic 
theories (some examples are discussed in section 6 of [44] ) , avoiding the conclusions of the 
Nielsen-Chadha theorem [45]. 

Said this, the behavior of the scalar correlator is quite surprising to us. Our naive 
expectation was that O would reproduce the propagator of a weakly coupled dilaton, that 
would lead to a different momentum dependence. Assume that, as naively expected, the 
low energy effective field theory description consists of the IR CFT plus an almost free 
dilaton. The dilaton does not couple directly to marginal operators, so the only allowed 
couplings to the CFT are to irrelevant operators, that would decouple at very low energy. 
The effective action at low energies in that case would be 

Cm ^ ClRCFT + ^dilaton + 9 ((C)) e A ' RT AlR . 

Here g{{0)) is the coupling of the irrelevant operator that drives the IR theory away from 
the fixed point in the RG flow. This coupling depends on the expectation value of the 
operator that breaks conformal invariance. Using that (OirOir) ~ {y / —q 2 ) 2A ' R ~ d , the 
dilaton propagator would take the form 



q 2 + (v / -g 2 ) 2Ajfl_d 

Clearly, this is different from the scalar correlator (OO) we have computed. In our case, the 
low-momentum behavior of the correlator coincides with that of an operator of dimension 
d - A IR in the IR CFT. 

(OO) ~ (V z <?) 2(d_A/ii)_d - 

This is the dimension of the coupling g((0)). We do not have a formal argument in field 
theory that explains this behavior, heuristically it looks like fluctuations of the expectation 
value lead to fluctuations of the coupling in the effective IR theory. Note that there is a 
strong IR divergence, the space-dependent Euclidean correlator grows like a power-law at 
long distances: 

(O(x)O{0)) ~ \x\ 2AlR . 

This is even stronger than the logarithmic divergence of a massless scalar field in two 
dimensions, that precludes the formation of a condensate [34]. Taken at face value, it 
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implies that there is no unitary field theory dual to the class of models we have chosen. This 
behaviour may be avoided in different models with no global definition of a superpotential 
or with more than one scalar field. Recall that having a globally defined superpotential is 
not a generic situtation. The power law we observe in the scalar correlator seems to be 
related to the coefficient of the quadratic term of the superpotential at the IR critical value 
4> m , which for a single scalar field in the bulk is necessarily negative. However, when there 
are more scalar fields classical trajectories in field space could connect two critical points 
of the superpotential along directions of positive curvature, a cartoon is given in figure |2[ 
Another possibility is that there is no global definition of the superpotential, so close to 
the IR critical value 4> m the local form of the superpotential has a positive coefficient in 
front of the quadratic term. 



Figure 2: Cartoon of a trajectory between two critical points of the superpotential in a two- 
dimensional space. The trajectory approaches the critical points along the directions of positive 
curvature. 

It would be interesting to study more general models with one or several scalar fields 
and see if the singular behaviour can be avoided. These models are also interesting for other 
reasons, for instance they can also be used to study the spontaneous breaking of Abelian 
or non-Abelian global symmetries. A particularly intriguing fact is that, to our knowledge, 
there are no examples in string theory models of geometries interpolating between two AdS 
spaces where the flow is dual to a spontaneously broken theory. The known models are 
either singular or involve explicit breaking. 

Summarizing, our result for the correlators leads to the following possibilities: 

• For some reason the low momentum expansion breaks down and there are additional 
terms in the scalar correlator that makes it less singular, maybe introducing a q 2 
term in the denominator. We checked corrections to the next few orders and could 
not figure out how this would happen. We also checked the method in the Coulomb 
branch solutions where the analytic solutions are known and found no problem there, 
but it would still be desirable to have an analytic example where the horizon geometry 
is AdS and the low momentum expansion is not necessary. 

• The singular behavior is an artifact of the large- N, strong coupling approximation 
that is implicit in the holographic approach and a proper treatment will avoid this 
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problems. This would imply that the supergravity approximation breaks down, which 
is difficult to understand from the bulk perspective, since the classical solution is 
completely smooth and the curvature is small everywhere. 

• The simplest gravitational models of a flow to a fixed point due to spontaneous 
breaking of symmetry have no unitary holographic duals. This could in principle be 
avoided in more general cases, even though there are no known examples. 

We would like to explore these issues in the future. 
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A. Equations of motion of fluctuations 

It is convenient to write the equations of motion ((O]) in terms of the extrinsic curvature 



^ = -^ = -8^, K* v =grK oa , = T$ v , K = gTK ia ,. (A.l) 
The Einstein equations take the form (primes denote radial derivatives) 

K> + K" p Kf> a = -Wf - JLv, (A.2) 

d^K - V a K% = -(j/dnh (A.3) 

iV + K'^ - 2K m K a u + KK^ = -dy/bdvb - J^j9n»V, (A.4) 

0" + K 0' + sTV^ - d^V = 0, (A.5) 

where hats denote covariant derivatives with respect to g^. It is convenient to choose a 
slightly different set of equations. For this, we will use the formulas 

K' = {sTK^y = gTK'^ - 2K«K au , (A.6) 

K'^ = (g^Kl)' = g^K a J + IK^K™ . (A.7) 

Now we take the trace of QA.4|) and substract ( |A.2j ) to get a new equation with only one 
radial derivative. Our last equation will be QA.4Q , but with one index raised 

K%K^ — K 2 + (0') 2 = 2V + R + gTd^h (A.8) 

d^K - V a K a ^ = -418^, (A.9) 

(K^Y + KK» V = -BP V - 8^8 u ct> - j^jKV, (A.10) 

(f/' + K^f + g^V^-d^V = 0, (A.ll) 
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In this case ( |A.£| ) and ( A. ICQ are constraints associated to the gauge fixing. We now expand 
the equations to linear order in fluctuations 



We will need the extrinsic curvature 



+ 



K% = A'5£ + ~<, K = dA' + lh', 



(A.12) 



(A.13) 



where primes denote radial derivatives, indices are raised with the flat metric and h 
h 1 ^ = ff^hyv. We will also use the Ricci tensor to linear order 

1 



-e 
2 



-2A 



Uh% - & t d x h\ - d v d x h" x + &*d v h 



(A.14) 



We now expand in Fourier modes, each mode admits the following decomposition in trans- 
verse and longitudinal parts, traceless and not traceless, as defined in (4.7) and ( |4.9| ). The 
equations of motion are 



= (d - l)A'ti + 2dVtp - 2(j)' (p' - q 2 e 



-2 A 



hn 



= qpti - q a h a J + 2q ft <// (p, 
Q = h tM J' + dA'h^ , + ,)': ( .1.'// + 



-dVp 



— e 



fq P V au 



d-l 

V qaTp t +»"» l, » a '' 



rf p )h. 



n // I j At t a2 T / 2 -2A , rOi/ 

= (p + dA ip — a Vip — q e ip + -^-h . 



(A.15) 
(A.16) 



(A.17) 
(A.18) 



Projecting (|A.16 ) with Pij? , gives the equation hT^ = 0, so we can set this mode to zero. 
Projecting the same equation with P£ v gives 



= h' T + 2<f/ <p. 



(A.19) 



Projecting with II Q ^ (A.17) gives a decoupled equation for the transverse traceless mode 



= + dA'h T j' - q 2 e~ 2A h T J. 



Projecting the same equation with either P T P or P L P gives the two equations 



= h' T + dA'h' T + (d- l)A'ti + 43V<p - q 2 e 



-2 A 



hn 



= h'[ + dA'h' L + A'ti + - — -dVip - q 2 e- 2A h T . 



(A.20) 

(A.21) 
(A.22) 



It turns out that not all the equations are independent. If we subtract ( [A. 15 ) from ( A.21] ) 
and use ( A. IE ) to remove the dependence on hx we find that the result vanish identically. 
Therefore, from (|4.12| ), ( |A.21[ ), ( A.19 ) and ( |A.22| ) we only need to keep three equations. 
We will combine flATp , flAgg ) and JO^ ) such that we find an equation that depends 
only on h and ip: 



= h" + 2A'h' + 



d-l 



(A.23) 
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We can determine (p and h from (|A.18 ) and ( |A,23 ), h' T from ( A.igj ) and h' L = h' — h' T . For 
the additional constraint ( A.15| ) to hold, it is enough to check that it is satisfied at some 
fixed value of the radial coordinate. Indeed, taking a derivative with respect to the radial 
coordinate and using the other equations, including the equations for the background, one 
can check that the result vanishes identically. 

Let us define now b! = e~ 2A H, we are left with the equations 



-2A 



H' + 



d-1 



dVip + 4$^', = ip" + dA'tp' - d 2 V<p - q 2 e- 2A y + —e~ 2A H. 



(A.24) 



B. The matching procedure 

In this section we explicitly show the existence of an overlapping region between the bound- 
ary and the near- horizon solutions. 



The equation for G ( 4.34 ) can be written in the following form 

x 2 d 2 G + a{x)xdG + p(x)Q 2 x$G = 



(B.l) 



where x = </>o — 4> m {4>m is the radial location of the horizon). The functions a(x) and /3(x) 
can be expanded around the horizon x = 



with 



a(x) = «o + otix + ct2X 2 + ... 
P{x) = /3 + /3ix + /3 2 x 2 + ... 

d-1 



A 

Lir\ 



and 



A) 

A = d - A IR < 



(B.2) 
(B.3) 



(B.4) 
(B.5) 

(B.6) 



Very close to the horizon we can approximate a = ao and (3 = With this approxi- 
mation the solution is (up to an overall constant that we drop) 



Gc, 



n f(i-«o) 
AQx^j #-A (1 _ ao) 



(B.7) 



where H is the Hankel function. In order to carry out the matching procedure we now 
have to expand this solution for Qx* <C 1. This expansion takes the form 

G = a 1 + ci 

+ b 1 + d x (jpoXQxt ) + ... ( y/PoXQx? 



7 /3oAQxa +... 



A(l-ao) 



(B.« 
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Our goal is to chek whether subleading contributions in the near-horizon expansion can 
give terms with the same dependence on x as the leading terms in the expansion of Gq. 
The first correction for the solution around the horizon is determined by the equation 



x 2 d 2 Gi + aoxdd + p Q Q 2 x^Gi 



aixdGo + (3iQ 2 x x Go 



(B. 



The exact solution for G\ will be some complicated superposition of Hankel functions. 
However, we are only interested in the expansion of the solution for large values of x, so 
we can solve this equation order by order. The solution takes the form 

#Q 2 xi + #Q 4 xi + .. 



ax 



+ bx 



/ I3q\Qx~: 



A(l-ao) 



(B.IO) 



where # represents numerical coefficients (with no dependence on the momenta Q). In a 
similar manner one can solve for higher orders and the result will take a similar form (with 
the coefficients in front of each of the series above replaced by ax n and bx 11 respectively, 
for the nth order). 

We conclude that the solution in the vicinity of the horizon takes the form 
G = a 1 + ci(x)Q 2 x^ + c 2 {x)Q A x^ + ... 



+ b 



do (a?) + di(x)Q 2 x* + d 2 (x)Q A x- + ... 



%\Qx* 



A(l-ao) 



(B.ll) 



where each of the functions above {cj,dj} can be expanded as a series in x and does not 
depend on the momenta Q. In particular, we find that do(0) = 1. 

By examining the power structure of the expansion above, one notices that some terms 
appear in both series. For instance, the order v term of the second series will mix with the 
order [v — of the first series. However, the leading terms in each series do not mix. 

9 — 

Therefore, when taking the limits i<1 and Q x* <C 1 only the near-horizon solution 
(Go) will survive and there would be no contributions from higher orders. That shows that 
there is a region where the boundary solution ( |4,36| ) and the near-horizon solution Q4.3£| ) 
overlap. 



C. Matching procedure for Coulomb branch solution 

In this section we check the matchig procedure at low momentum for the Coulomb branch 
solutions. Our conventions are different from those in [25], both for the gauge fixing and 
the definition of the scalar action. In our conventions, the superpotential is 

v + 2 



W 



,1/3 ' 



,V3<f> 



(C.l) 



C.l Tensor fluctuation 



The equation of motion for the tensor fluctuation is the same as (4.11), using v as variable 
we can also write it as 



v(v - l)d 2 H T + (v- 2)d v H 7 



a 2 

5 ttT 



Av 



-H 1 



0. 



(C.2) 
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In this coordinate, v = is the horizon (which is singular) and v = 1 the boundary. The 
exact solutions are 



H T = Cov~ a 2 F 1 



a a 

--,--,2-a; V 



+ C 1 v a ~ 2 2 F 1 



a a 

2" 1 '2" 1 ' a;W 



(C.3) 



where a = \/l + Q 2 + 1. Imposing regularity at the horizon sets Co = and C\ = Hq, 
The expansion of the regular solution close to the boundary at low momentum gives 



(C.4) 



So A\ can be identified as the source of the energy-momentum tensor up to kinematical 
factors. 

Let us now see how this behavior is captured by the matching procedure. First, we 
take the near-horizon limit v — > 0, where the equation of motion (|C.2|) takes the form 



2 a 2 
8 2 H T + -8 V H T - -^H T = 0. 



The solutions to this equation are 



H 1 = A v2 1 +A x v 

Expanding at low momentum Q 2 log «<1 gives 

A 



H T ~ -± [1 + Q(Q 2 log v)] + A [1 + 0(Q 2 log «)] . 



(C.5) 



(C.6) 



(C.7) 



We now consider the q 2 — > limit, where the solution has the form ( |4.16| ), and expand 
for small v 



1 



Ho 1 - -Hi y- + \ gv )+0(Q 2 ) 



Including the correction to next order we have 



1 



Q 



H 1 * H6 + ( —Hi - 2-{H? - Hi) ) f- - -Hi logv ( 1 + X7 



2 1 4 
We can match the two expressions if 



1 1 



v 2 



Q' 2 



4l> 



A = Hq,Ax = --Hf - -%-(Hi - Hq), 



(C.8) 



(C.9) 



(CIO) 



and Hj oc Q 2 so the expansion is well defined. The regularity condition A\ = implies 



H 



o - 



(C.ll) 
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to leading order in the expansion. Therefore, the solution becomes 



Q 2 



H 1 ~ H< [ 1 - ^-logv + 0(Q* 



^ HZ I 1 



- 1) + f> - 1) 2 + o(q 4 ; 



(C.12) 



in agreement with the expansion of the exact solution (C.4). 
C.2 Scalar fluctuation 



The equation of motion for the scalar fluctuation is given by (4.31). In terms of the v 
coordinate it is given by 



3 13-lOv 2 3q 2 (v-l)+8(v 2 -2v+4) 
& ^ + 3v _ 3v 2 9 vV + Ld vf 



X2(v-l) 2 v 2 

q 2 (-2v 2 +v+l)-24v 
^ 12(v-l) 3 « a ^ 



0. 



(C.13) 



Close to the boundary u — >• 1 there is a regular singular point, so this equation has three 
independent solutions 



tp ~ a (v -l) + ai(v- 1) log(l -v)+ (a 2 - yj (v - l) 2 



(C.14) 



The solution proportional to ao has the asymptotic behaviour of the gauge solution (p = 
dW, that is an exact solution to the equations of motion. The other two solutions can be 
obtained from the computation of the gauge-invariant combination ( |4.33| ) , that obeys the 
equation (4.34). In terms of the v coordinate, it becomes 

„2 



d 2 v G+^^d v G + 



v + 2 



4(u - l)v 2 



G = 0. 



The exact solutions are 
G = A v 1 ~ a / 2 2 F 1 



a a 



-, 2 — a; v 



+ A x v a ' 2 2 F l 



a a 
2" 1 '2' a;W 



(C.15) 



(C.16) 



where a = \fl + q 2 + 1. Imposing regularity at the horizon sets Aq = 0. For the regular 
solution, the expansion close to the boundary at low momentum is, to leading order 



^-^\ + 0{q 2 ) + (v-l) 



I + 0(q 2 ) - ^ log(l -v) + ^(v- I) 2 log(l - v) 



. (C.17) 



Let us know show how the matching works. First, we approximate ( |C.15|) by the 
near-horizon form 



d 2 G - -j-yG = 0. 

4 V 2 



(C.18) 



Solutions to this equation are 



G~ Aqv 1 -^ 2 + A lV ^ 2 , 



(C.19) 
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the regularity condition being Aq = 0. Expanding for low momentum gives 

G~^-v(l + £]ogvY (C.20) 
We now solve the equation in the low momentum limit 

G ~ Co + -^cj 2 + 3(V " 1} log(1 " V) + (2 " 3v) bg V 



v + 2 4 V 4(w + 2) 

<?V - l)log(l -v) - ulogv + logw + I s . 

T Cl I —2 I • (c ' 21) 



Expanding close to the horizon, one gets 



2 J \ 4 i 
+ ( - + 4C ) + ^ Iom rCk ', + SC '„ ) j . (C.22) 



Setting C = -C\/2 and Ci = -2A X matches with ( p.2CQ 



Expanding now ( C.21| ) close to the boundary, 



G~ A x [\ + {v-l)(l-f-\og{l-v) + ^-(v-l) 2 \og{l-v) ) ) . (C.23) 



3 v 'V9 12 36 



This is in agreement with ( p,17| ) , therefore the matching procedure fixes the leading terms 
in the low momentum expansion. 

D. RG flow toy models 



J m)i 



In order to have an AdS metric when <f> ~~ 0m , it is necessary that W'((j)) ~ (4> — <ft r . 
so the superpotential has a critical point W'(<f> m ) = 0. 8 Then, in order to have a soliton 
solution between <p = and (j) = 4> m , we will require that the superpotential satisfies the 
following conditions (for (f) m > and 5 = Lyy/LiR > 1) 

W(0) = ^-—^, W'(0)=0, W"(0) = ^>0, (D.l) 
W(<p ) = ^, W'(<Po) = 0, W"(fo) = ^Z^m <0 . (D.2) 

J->IR L IR 

There are many such possible superpotentials, so we need some additional criteria to make 
our choice. In supergravity theories this kind of superpotentials can be found after doing a 
dimensional reduction in some internal manifold (the scalars correspond to the moduli) and 



8 Notice that critical points of the superpotential are always critical points of the potential, since V'(<f>) = 

w'(4>) (w"(4>)-^xW(ct>)). 
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are usually a sum of exponentials. Inspired by this fact, we will make our superpotential 
also a sum of exponentials. In order to find two critical points we need at least three 
exponential terms, a possible "minimal" superpotential is then 



W{4>) = W e a4> + 



+ k 2 



2a<t> 



1 + k + k 2 2(1 + k + k 2 ) 



(D.3) 



By construction, W((f)) has two critical points, one at <j) = and the other at 
— l/alog(fc). 

Using the scalar field as radial coordinate, 



ds 2 

the coefficients of the metric are 
1 1 



+ e 



2A(<f>) 



r\ ilv dx iL dx h 



(W'(4>)f a 2 W 2 



(k + k 2 ) 2a( p 



and 

A{4>) = A + 



1 



a(d-l) 2a 2 (d-l 
log ( 1 — ke 



l + k + k 2 2(1 + k + k 2 ) 



4 + k + k 2 



(D.4) 



(D.5) 



atj> 



2-k-k 2 

4 + 7k + 4k 2 



log ( e aq> - 1 

2 

log ( 1 + (fc + l)e a? 



;2 + A;)(l + 2it) 



(1 - k)(2k + \ 

where Aq is an arbitrary integration constant. The radial coordinate is 

(1 + k + k 2 ) r, a k(k + l 



(D.6) 



r(<P) 



W 
1-k 2 
(l + 2k) 



a 2 W (2 -k-k 2 ) 
log ( 1 + (k + l)e a < 



log e 



1 



(2k + 1) 



log (1 - ke^ 



(D.7) 



where rg is the integration constant. 



One can satisfy all the conditions (D.l) if 
d-ll + k + k 2 ~ 



Close to 



0. 



L uv 4 + k + k 2 '' 



r ~ — log((p) 



4 + k + k 2 



k - k 2 V 2(d - 1) ' 



< k < 1. 



1 



A 



log( 



~ e - A uvr/L uv ^ r _^ ^ 



e 2A ^ e 2r/L uv _ 



(D.8) 

(D.9) 
(D.10) 



This shows that the metric indeed becomes AdS. Furthermore, the profile of the scalar 
field corresponds to the normalizable solution, which means that in the dual field theory 
the scalar operator has a vacuum expectation value but there is no explicit breaking of 
conformal invar iance. 
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The ratio between the UV and IR AdS radii is 



L IR _l_ fc(4 + fc + fc 2 ) 

L uv 5 l + k + Ak 2 ' K ' ' 

The IR dimension is 

A J r w ,/ r< x \ _ J i A (4 + fc + fc 2 )(l + fc-2fc 2 ) 

We have constructed a one-family parameter of soliton solutions for any given d and d > 
A uv > 0. 
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